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PREFACE 


This  investigation  was  sponsored  by  the  Defense  Nuclear  Agency 
under  Subtask  Y99QAXSB211,  "Penetration,"  Work  Unit  l8,  "Penetrator 
Summary."  This  study  was  conducted  by  personnel  of  the  Soils  and  Pave- 
ments Laboratory  (S&PL),  U.  S.  Army  Engineer  Waterways  Experiment  Sta- 
tion (WES),  during  January  through  December  1977,  under  the  general 
supervision  of  Messrs.  J.  P.  Sale  and  R.  G.  Ahlvin,  Chief  and  Assistant 
Chief,  S&PL,  respectively.  Mr.  R.  S.  Bernard  formulated  the  theory 
and  Mr.  D.  C.  Creighton  implemented  the  computer  analysis,  both  under 
the  direct  supervision  of  Dr.  J.  G.  Jackson,  Chief  of  the  Soil  Dynamics 
Division  (SDD),  S&PL.  Drs.  P.  F.  Hadala  and  B.  Rohani,  SDD,  provided 
technical  guidance.  Mr.  Bernard  and  Mr.  Creighton  prepared  this  report. 

COL  J.  L.  Cannon,  CE,  was  Commander  and  Director  of  WES  during 
the  period  of  research  and  report  preparation.  Mr.  F.  R.  Brown  was 
Technical  Director. 
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CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Si) 
UNITS  OF  MEASUREMENT 


U.  S.  customary  units  of  measurement  used  in  this 
verted  to  metric  (SI)  units  as  follows: 


Multiply 


feet  per  second 

inches 

pounds  (force) 
per  square  inch 

inch-pounds 

pounds  (mass) 

pounds  (mass) 
per  cubic  foot 

pounds-square 

inches 

slugs  per  cubic 
foot 

square  inches 


0.301+8 

0.025!+ 
689I+.  757 

0.1130 

0.1+535921+ 

16.0181+6 

0.0002926 

515.50336 

0.00061+516 


report  can  be  con- 

To  Obtain 

metres  per  second 

metres 

pascals 

newton-metres 

kilograms 

kilograms  per 
cubic  metre 

kilograms-square 

metres 

kilograms  per 
cubic  metre 

square  metres 
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NON-NORMAL  IMPACT  AND  PENETRATION:  ANALYSIS  FOR 


HARD  TARGETS  AND  SMALL  ANGLES  OF  ATTACK 


CHAPTER  1 
INTRODUCTION 


1.1  BACKGROUND 

There  now  exist  several  techniques  that  are  useful  for  making  pene- 
tration predictions,  provided  the  projectile  survives  the  impact  and  re- 
mains stable  thereafter.  The  projectile  orientation  and  motion  are  usu- 
ally assumed  normal  to  the  target  surface,  but  perfectly  normal  impact 
seldom  occurs  in  practice.  A small  angle  of  attack  is  ordinarily  to  be 
expected. 

The  Defense  Nuclear  Agency  (DNA)  recently  initiated  an  investiga- 
tion of  non-normal  impact  and  penetration,  with  emphasis  on  small  attack 
angles  (less  than  10  degrees1).  As  a part  of  this  effort,  DNA  asked  the 

U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES)  to  develop  a non- 

2 

normal  penetration  analysis  based  on  the  Cavity  Expansion  Theory  (CET). 
AVCO  Corporation  (Reference  l)  had  previously  developed  a computer  code 
for  analyzing  projectile  motion  in  three  dimensions,  using  an  empirical 
Differential-Area  Force  Law  (DAFL)  to  specify  the  surface  loads  on  the 
projectile.  The  AVCO  version  of  the  DAFL  code  input  contains  a set  of 
user-selected  empirical  parameters  describing  the  target,  and  it  was 
felt  that  the  CET  analysis  might  offer  a simpler,  more  rational  target 
description. 

1.2  PURPOSE  AND  SCOPE 

The  purpose  of  this  study  was  to  extend  the  WES  CET-based  penetra- 
tion theory  for  normal  impact  (Reference  2)  to  include  non-normal  impact. 


A table  of  factors  for  converting  U.  S.  customary  units  of  measurement 
to  metric  (SI)  units  is  presented  on  page  6. 

The  Cavity  Expansion  Theory  is  discussed  in  Appendix  A. 
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With  the  CET  analysis  in  place  of  AVCO's  force  distribution,  a modified 
version  of  the  DAFL  code  was  used  to  make  calculations  for  comparison 
with  test  data.  Sandia  Laboratories  and  AVCO  Corporation  have  each  con- 
ducted a series  of  instrumented  reverse  ballistic  tests  (RBT's)  for  DNA 
at  attack  angles  between  0 and  10  degrees.  Comparison  with  these  tests 
shows  the  CET  analysis  to  be  tentatively  acceptable  as  a rigid-body  pre- 
diction technique  for  hard  targets  and  small  angles  of  attack. 
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CHAPTER  2 


THEORY 

2.1  PROJECTILE  MOTION  IN  TWO 
DIMENSIONS 

Previous  CET -based  penetration  analyses  (Reference  2)  have  ad- 
dressed only  the  axial  motion  of  the  projectile,  obtaining  expressions 
for  the  total  resisting  force.  The  non-normal  problem  in  two  dimen- 
sions, however,  involves  three  degrees  of  freedom:  one  rotational, 
and  two  translational  (Figure  2.1). 

Consider  a rigid,  deeply  buried  projectile  (no  free-surface  ef- 
fects), whose  motion  lies  in  the  XZ-plane.^  The  rotation  is  described 

• • 
by  the  angular  velocity  0 ; the  translation  is  described  by  X and 

Z , which  are  the  X-  and  Z-velocity  components  of  the  center  of  gravity 

(CG).  ’ For  analysis,  it  is  convenient  to  express  the  CG  velocity  in 

terms  of  its  lateral  and  axial  (x-  and  z-)  components  V and  V , 

x z 

respectively: 


V = X cos 

0 

• 

- Z sin 

0 

(2.1) 

X 

V = X sin 

0 

+ Z cos 

0 

(2.2) 

For  purely  axial  motion  (V  =0,0=0),  the  force  distribution  on  the 
projectile  is  symmetric,  and  the  net  force  is  in  the  z-direction.  The 
introduction  of  rotation  (0^0)  and  lateral  motion  (V  / 0)  destroys 
this  symmetry,  and  simplifying  assumptions  are  necessary  in  order  to 
calculate  the  resulting  (asymmetric)  force  distribution. 

Whenever  the  lateral  and  rotational  velocities  are  small 

1 

2 

3 
k 


The  XZ-coordinate  system  is  fixed  in  the  target  (Figure  2.1). 

A dot  above  any  quantity  indicates  differentiation  with  respect  to 
time. 

For  convenience,  symbols  and  abbreviations  used  in  this  report  are 
listed  and  defined  in  the  Notation,  Appendix  D. 

The  xz-coordinate  system  is  fixed  on  the  projectile  axis,  with  its 
origin  at  the  CG  (Figure  2.1). 


9 


(V  <<  Vz  , L6  <<  Vz  where  L is  projectile  length),  their  contribu- 
tions to  the  force  distribution  are  assumed  to  be  independent  of  the 
axial  motion.  Thus,  defining  the  angle  of  attack  as 


the  force  contributions  due  to  axial,  lateral,  and  rotational  motion 

can  be  defined  separately  when  tan  a <<  1 and  6 <<  V /L  . ^ 

z 

Initially  it  is  assumed  that  the  entire  projectile  is  in  contact 
with  the  target.  This  assumption  will  later  be  modified  when  the  ef- 
fect of  wake  separation  is  investigated.  Whether  or  not  there  is  con- 
tact at  a given  point,  the  following  assumption  will  remain  in  effect 
throughout  the  analysis:  When  any  point  on  the  projectile  surface 
moves  away  from  the  adjacent  target  material,  the  stress  at  that  point 
is  zero. 

2.2  FORCE  DISTRIBUTION  DUE  TO 
AXIAL  MOTION 


For  purely  axial  motion,  the  axial  resisting  force  can  be  calcu- 
lated by  assuming  that  only  normal  stresses  act  on  the  projectile  sur- 
face.^ Assuming  further  that  the  axial  force  is  unchanged  for  small 

attack  angles,  the  following  equation  is  used  to  approximate  the  com- 

7 

pressive  normal  stress  due  to  axial  motion  only: 

o=a+  V «/po  ( 1 + sin  n)  (2.1*) 

» S Z yf  S 

where 

5 • 

The  0 criterion  will  always  be  satisfied  for  stable  projectiles 
g when  tan  a <<  1 . 

The  analysis  in  Reference  2 contains  both  normal  and  tangential 
stresses.  For  hard  targets  such  as  rock  and  concrete,  however,  the 
j calculated  tangential  stresses  are  small  enough  to  be  ignored. 
Equation  2.1*  is  an  assumed  stress  distribution  that  produces  an  ap- 
proximation for  the  net  axial  force  when  integrated  over  the  pro- 
jectile surface. 
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r 


= radial  compressive  stress  at  the  surface  of  a slowly  ex- 
panding spherical  cavity  in  the  target  material  (Appendix  A) 

p = target  mass  density 


tan  n 


slope  of  projectile  surface  at  a given  point,  with  respect 
to  axis  of  symmetry^  (Figure  2.1) 


The  x-  and  z-components  of 


a 

a 


are,  respectively. 


o'  = a cos  n cos  ip 
x a 


(2.5) 


o'  = o sin  n (2.6) 

Z Q, 

where  ip  is  the  azimuthal  angle  (Figure  2.2).  For  purely  axial  motion, 
with  complete  contact  between  target  and  projectile,  only  Equation  2.6 
contributes  to  the  net  resisting  force  (no  net  lateral  force).  Other- 
wise there  i£  a net  lateral  contribution  from  Equation  2.5. 

2.3  FORCE  DISTRIBUTION  DUE  TO 
LATERAL  MOTION 


For  small  attack  angles,  the  resisting  force  due  to  lateral  motion 

is  presumed  to  be  unaffected  by  the  axial  motion.  In  this  sense,  the 

lateral  motion  can  be  regarded  as  a "separate"  penetration  problem,  and 

it  is  assumed  that  the  lateral  motion  produces  only  lateral  resistance. 

Dropping  the  nose  shape  factor  (l  + sin  n)  and  replacing  V with 

z 

V , the  expressions  for  the  x-  and  z-components  of  the  compressive 

X 9 

stress  due  to  lateral  motion  are  assumed  to  be 

/ 


a"  = 
x 


+ V 

x 

- V 

x 


pa  ) cos  n cos  ip 


COS  n cos  ip 


V > 0 
x 

V < 0 
x 


(2.7) 


The  quantity  1 + sin  n , which  appears  in  Equation  2.U,  is  an  empiri- 
cal "nose  shape"  factor,  introduced  to  produce  a gradual  increase  in 
q axial  resistance  with  increasing  nose  bluntness  (sin  n). 
y Equations  2.7  and  2.8  are  analogous  to  Equations  2.5  and  2.6.  The 
nose  shape  factor  is  eliminated  in  the  lateral  stress  because  all  pro- 
jectiles of  interest  have  circular  cross  sections. 
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o' ' 
z 


0 


(2.8) 


Wherever  cos  4>  < 0 , the  lateral  translation  of  that  point  is  away 
from  the  adjacent  target  material,  and  so  a”  = 0 for  that  point  on 
the  projectile. 

2.b  FORCE  DISTRIBUTION  DUE  TO 
ROTATION 

Assuming  that  the  lateral  velocity  due  to  pure  rotation  z0  merely 
adds  an  additional  velocity  term  to  Equation  2.7,  the  x-  and  z- 
components  of  the  compressive  stress  due  to  rotation  are 

o^' ' = z0  y pog  cos  n cos  \l>  (2.9) 

0"  ' = 0 (2.10) 

z 

Wherever  z0  cos  \p  < 0 , the  rotation  is  away  from  the  adjacent  target 
material,  and  so  o^' ' = 0 for  that  point  on  the  projectile. 

2.5  COMPOSITE  FORCE  DISTRIBUTION 

Equations  2.5  through  2.10  can  now  be  combined  to  obtain  the  com- 
posite axial  and  lateral  force  distributions  on  the  projectile.  The 
net  axial  stress  is 


o =a,+a',+a,,,=o,  (2.11) 

z z z z z 

The  net  lateral  stress  is 

0 = o'  + 0"  + 0'”  (2.12) 

X X X X 

The  lateral  and  axial  force  components  (Figure  2.2)  acting  on  any  sur- 
face element  dA  are,  respectively, 

dF  = -0  dA  (2.13) 

x x 
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(2.11*) 


dF  = -o  dA 
z z 

where  dA  is  the  differential  surface  area.  The  lateral  and  axial 
force  components  acting  on  the  entire  projectile  are  then 

F = - / a dA 
x J x 

F = - / o dA 
z J z 

A 

where  A represents  the  entire  projectile  surface.  The  translational 
equations  of  motion  are 

MX  = F cos  0 + F sin  0 (2.17) 

x z 

MZ  = - F sin  6 + F cos  0 (2.l8) 

x z 

where  M is  the  projectile  mass  and  X and  Z are  the  components  of 
acceleration  in  the  X-  and  Z-directions , respectively.  The  projectile 
is  assumed  rigid,  and  the  rotational  equation  of  motion  is 

10  = / z dF  - / x dF  (2.19) 

A x JA  z 

where  I__  is  the  mass  moment  of  inertia  about  the  CG,  and  6 is 
the  angular  acceleration.  The  right-hand  side  of  Equation  2.19  is  the 
total  moment  exerted  on  the  projectile  about  its  CG. 

2.6  WAKE  SEPARATION  AND 
REATTACHMENT 

It  is  likely  that  some  parts  of  the  projectile  lose  contact  with 
the  target  during  penetration.  Part  of  the  nose  must  remain  in  contact, 
but  separation  can  occur  anywhere,  especially  on  the  aftbody.  While  a 
precise  analysis  of  the  separation  mechanism  itself  is  unmanageable,  it 
is  still  possible  to  assess  the  gross  effect  of  a specified  amount  of 
separation  upon  the  projectile. 

i 


(2.15) 

(2.16) 
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Consider  a projectile  with  no  rotational  or  lateral  velocity 
(5=0  , V = 0).  Separation  is  thought  to  occur  near  the  base  of  the 
nose  (Figure  2.3),  with  reattachment  somewhere  on  the  aftbody  (if  any- 
where). The  separation  and  reattachment  points  probably  change  with 
velocity,  and  the  width  of  the  separation  gap  may  vary  with  position. 
This  behavior  is  simplified  and  idealized  (Figure  2.M  in  the  follow- 
ing assumptions: 

1.  The  entire  nose  is  in  contact  with  the  target. 

2.  The  width  of  the  separation  gap  6 between  aftbody  and 
target  is  directly  proportional  to  the  local  aftbody 
radius  r . 

Assumption  1 is  reasonable  for  long  projectiles,  for  which  the 
nose  represents  only  a fraction  of  the  total  surface  area.  According 
to  Assumption  2,  a cavity  exists  in  the  target  aft  of  the  nose  (Fig- 
ure 2.1+) , such  that  6/r  remains  constant.  These  assumptions  have  no 
effect  whenever  V =0  and  0=0.  However,  if  V / 0 or  0^0, 
then  there  is  a distribution  of  local  lateral  velocity 

v = V + z0  (2.20) 

xx 

and  part  of  the  aftbody  surface  may  collide  with  the  cavity  wall  (which 
is  regarded  as  fixed).  Reattachment  is  defined  as  any  collision  be- 
tween the  aftbody  and  the  cavity  wall. 

The  problem  is  still  complicated  because  the  cavity  wall  will  not 
be  straight  if  / 0 or  0 / 0 . Instead,  it  will  be  a curved  cy- 
lindrical surface,  mapping  out  the  trajectory  of  the  nose  base.  Addi- 
tional simplification  is  necessary  if  the  analysis  is  to  be  tractable. 

Consider  the  two  views  of  the  idealized  wake  separation  shown  in 
Figure  2.h.  Suppose  that  the  cavity  is. almost  straight  between  points 
1 and  2,  but  that  point  1 (on  the  projectile)  is  moving  toward  the 
cavity  wall.  The  simplified  reattachment  criterion  is  as  follows:  if 
point  1 can  collide  with  the  cavity  wall  before  going  past  point  2,  then 
point  1 is  assumed  to  be  already  in  contact  with  the  cavity  wall.  Thus, 
the  condition  required  for  contact  is 


(2.21) 


where  s is  the  axial  distance  between  the  point  in  question  and,  the 
base  of  the  nose.  At  arbitrary  points  on  the  aftbody  (such  as  point  3), 
one  of  the  following  conditions  is  required  for  contact: 


V2  2 2 

r cos  + 2r6  t { - r 


cos  i|/  _ s w IT 

< — . ~ 2 * K 2 (2-22) 


v; 


2 2 2 

cos  <h  + 2r 6 + 6 + r cos  1 1>  s 

Vz  ’ 


| < ♦ < y (2-23) 


Wherever  this  condition  is  satisfied  (on  the  aftbody),  the  stress  is 
defined  by  Equations  2.11  and  2.12;  otherwise,  the  stress  is  set  equal 
to  zero  (no  contact). 

2.7  INTERNAL  FORCES  AND  MOMENTS 
IN  A RIGID  PROJECTILE 


The  external  forces  and  moments  create  internal  forces  and  moments 
within  the  rigid  projectile  itself.  According  to  D'Alambert's  princi- 
ple, the  x-  and  z-components  of  the  internal  force  at  any  point  z on 
the  axis  of  symmetry  are,  respectively. 


fx(z)  = | [dFjO  - vx  dM(£)] 


fz(z)  = | [dFz(0  - vz  dM(0] 


{2.2k) 


(2.25) 


where  z is  the  z-coordinate  of  the  nose  tip,  the  differential  func- 
n 

tion  dM(z)  represents  the  point  distribution  of  mass  along  the 
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projectile  axis  (x  = 0),  and  £ is  a variable  of  integration.  The 
differential  force-components,  dF^(z)  and  dF^z)  , are  given  by 
Equations  2.13  and  2.lh.  The  local  lateral  and  axial  accelerations 
are  given  by 


\ 

i 

i 


V 

X 


d 

dt  Vx 


7-  V + z6 
dt  x 


(2.26) 


v 

z 


d 

dt  Vz 


— V 
dt  z 


The  internal  (bending)  moment  at  z is 


(2.27) 


£=z  x=+r 

r(z)  = / (zn  - 5)  dl^U)  + / x dFz(x,z) 

C=zn 

[(z  - 5)  v - xv  ] dM(x,£) 

n x z 

£=z  x=-r 
n 


(2.28) 


The  differential  functions  dM(x,z)  and  dF  (x,z)  represent  the  dis- 

z 

tributions  of  mass  and  external  axial  force,  respectively,  throughout 
the  projectile. 

2.8  MODIFIED  AVCO  COMPUTER  CODE 

The  equations  for  the  x-  and  z-components  of  stress  (Equations  2.11 
and  2.12)  have  been  installed  in  a modified  version  of  the  AVCO  DAFL 
code,  replacing  the  AVCO  definitions  of  the  same  quantities.  Algo- 
rithms have  also  been  added  for  calculating  internal  force  (Equa- 
tions 2.2b  and  2.25)  and  bending  moment  (Equation  2.28)  in  a rigid  pro- 
jectile. Figure  2.5  shows  a block  diagram  for  calculations  using  the 
CET  analysis  and  the  modified  DAFL  code. 
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2.9  RESTRICTION  TO  HARD  TARGETS  AND 
SMALL  ANGLES  OF  ATTACK 

The  theory  in  Sections  2.1  through  2.7  is  intended  only  for  hard 
targets  such  as  rock  and  concrete.  Furthermore,  since  the  analysis 
uses  purely  axial  motion  as  a point  of  departure  from  normal  impact, 
calculations  should  he  restricted  to  small  angles  of  attack  (a  10  de- 
grees), for  which  the  projectile  motion  is  predominantly  axial. 


Figure  2.1  Projectile  motion  in  two  dimensions 


Figure  2.2  Three-dimensional  viev  of  projectile. 


Figure  2.3  Wake  separation  and  reattachaent  for 
purely  axial  motion. 


Figure  2.5  Block  diagram  for  non-normal  impact-and-penetration 
calculations  using  CET  analysis  and  AVCO  DAFL  code. 
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CHAPTER  3 


CALCULATIONS  AND  TEST  DATA 

3.1  AVAILABILITY  OF  NON-NORMAL  TEST 

DATA 

Direct  measurements  of  the  lateral  forces  incurred  by  a projectile 
during  non-normal  impact  and  penetration  are  scarce.  In  fact,  with  the 
exception  of  tests  recently  conducted  by  Sandia  Laboratories  (References 
3 and  U)  and  AVCO  Corporation^"  (References  5 and  6),  lateral  force- 
related  data  are  virtually  nonexistent.  These  tests  seem  to  offer  the 
only  benchmarks  for  checking  the  theory  formulated  in  Chapter  2,  but  the 
available  data  are  too  few  to  establish  the  general  credibility  of  the 
theory.  More  tests  are  needed  under  different  impact  conditions,  and 
the  test  results  that  do  exist  are  sufficient  only  for  tentative 
verification. 

3.2  LIMITATIONS  AND  SIMPLIFICATIONS 

IN  THE  CET  ANALYSIS 

The  CET  analysis  is  intended  for  calculations  in  which  the  initial 
angle  of  attack  is  less  than  10  degrees.  The  primary  assumption  is  that, 
for  small  attack  angles,  the  lateral  resisting  force  can  be  found  by 
linear  superposition  of  axial,  lateral,  and  rotational  effects. 

It  will  be  shown  that  aftbody  wake  separation  has  little  effect  in 
the  calculations  if  the  assumed  separation  gap  width  6 is  small  com- 
pared to  the  local  aftbody  radius  r . For  hard  targets  it  is  likely 
that  the  projectile-target  separation  is  small;  thus  the  comparison  of 
calculations  and  test  data  will  be  made  on  the  assumption  that  6/r  ~ 0. 
The  parametric  effect  of  varying  6/r  in  the  calculations  will  be  shown 
separately. 

The  quantity  ag  , which  appears  in  Chapter  2,  is  discussed  in 


Personal  communication  dated  10  May  1977  from  P.  J.  Grady,  AVCO  Corpo- 
ration, Wilmington,  Mass.,  to  MAJ  D.  Spangler,  Defense  Nuclear  Agency, 
Washington,  D.  C. 
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Appendix  A.  Defined  in  terms  of  the  target  shear  strength,  links 

the  CET  to  the  penetration  theory.  For  calculations  in  hard  materials 
(Reference  2),  it  has  heen  found  acceptable  to  assume  that  the  target 
obeys  the  von  Mises  failure  condition.  The  strength  parameters  (Appen- 
dix A)  then  become 


(P1  = 41 2 = 0 


(3.1) 


C1  = C2  " 2 Y 


(3.2) 


where  Y is  the  unconfined  compressive  strength.  The  expression  for 


a reduces  to 

s 


o = f Y(1  + In  I.) 
s 3 1 


(3.3) 


where  1^  is  the  rigidity  index,  given  by 


_ _ 2E 
1 3Y 


(3. U) 


and  E is  Young's  modulus  for  the  target. 


3.3  SANDIA  REVERSE  BALLISTIC  TESTS 


3.3.1  Test  Description 

Sandia  Laboratories  has  conducted  a series  of  RBT's  for  both  normal 

2 

and  non-normal  impact.  In  each  of  the  experiments,  a sandstone  target 
was  mounted  on  a rocket  sled  and  accelerated  to  a velocity  of  about 
1500  ft/s  before  striking  a stationary  1.9-inch-diameter  projectile. 

Two  of  the  tests  were  conducted  at  attack  angles  of  3 degrees;  the  im- 
pact conditions,  projectile  parameters,  and  sandstone  properties  for 
these  tests  are  shown  in  Figure  3.1. 

The  projectiles  were  instrumented  with  strain  gages  and  with  lateral 


The  targets  were  hexagonal  prisms,  1*5  inches  across  and  30  inches  deep 
(in  the  direction  of  motion). 
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ar.d  axial  accelerometers  (Figure  3*2),  and  data  were  recorded  for  the 
first  body  length  of  penetration.  The  outputs  of  the  accelerometers  (at 
stations  F and  R)  and  the  strain  gage  (at  station  S)J  can  be  compared 
directly  with  the  predictions  of  the  CET  analysis. 

3.3.2  Accelerometer  Data 

Figures  3.3  and  3.*+  show  comparisons  of  the  calculated  and  measured 
acceleration  records  for  the  forward  accelerometer  location  (station  F). 
The  axial  deceleration  calculation  (Figure  3.3)  reproduces  the  experi- 
mental results  about  as  accurately  as  can  be  expected  for  a rigid-body 
calculation.  The  lateral  acceleration  records  for  this  station  (Fig- 
ure 3.M  are  more  difficult  to  appraise.  The  CET  calculation  reproduces 
the  maximum  acceleration  levels,  but  not  the  instantaneous  details. 

Figures  3-5  and  3.6  show  the  calculated  and  measured  acceleration 
records  for  the  rear  accelerometers  (station  R).  Again,  the  axial  de- 
celeration calculation  is  in  reasonable  agreement  with  the  experimental 
results  (Figure  3-5).  Figure  3.6,  however,  offers  a more  dramatic  com- 
parison of  lateral  acceleration  records  than  does  Figure  3.b.  The  ex- 
perimental acceleration  records  show  about  the  same  magnitude  and  dura- 
tion as  the  calculated  results,  but  the  peaks  occur  later  than  the 
calculated  results.  This  delay  is  apparently  caused  by  nonrigid  projec- 
tile response  to  the  rapidly  varying  applied  loads  (Appendix  B). 

3.3.3  Strain  Gage  Data 

In  the  CET  analysis  it  is  assumed  that  the  projectile  is  perfectly 
rigid,  but  actually  strains  are  produced  in  the  projectile  by  the  inter- 
nal force  f and  bending  moment  T . In  the  plane  of  motion,  the  ex- 
z 

pression  for  the  maximum  axial  strain  on  either  side  of  the  projec- 

tile is 


ea  “ E A 
P s 


V 

E I 
p s 


(3.5) 


where 

E^  = Young's  modulus  for  projectile 


There  were  strain  gages  at  several  locations,  but  the  one  nearest  the 
CG  has  been  chosen  for  discussion. 
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= area  of  projectile  cross  section  S 

rg  = projectile  radius  at  section  S 

I = area  moment  of  inertia  of  projectile  cross  section  S 

If  e is  small  and  f and  r are  slowly  varying  functions  of  time, 

El  2 

then  Equation  3.5  is  an  exact  formula  for  the  strain  gage  output. 

Figure  3.7  shows  a comparison  of  the  calculated  and  measured  axial 

strains  for  strain  gage  S (Figure  3-1).  The  calculated  and  experimental 
strain  levels  are  about  the  same,  but  the  experimental  results  are  de- 
layed, as  were  the  lateral  acceleration  data.  Again,  the  delay  seems  to 
be  caused  by  the  rapidly  varying  applied  loads  (Appendix  B). 

3.U  PARAMETER  STUDY 

The  real  test  of  the  CET  analysis  lies  not  in  the  prediction  of  in- 
dividual test  results,  but  in  the  prediction  of  the  results  (or  trends) 
observed  in  many  tests.  The  data  required  for  such  a verification  are 
as  yet  unavailable,  but  it  is  still  worthwhile  to  investigate  the  im- 
portance of  different  parameters  in  the  calculations  themselves.  Accord- 
ingly, a brief  study  has  been  made  in  which  the  attack  angle,  the  impact 
velocity,  the  target  strength,  and  the  wake  separation  were  varied.  The 
Sandia  test  conditions  were  used  as  a baseline,  and  the  projectile  param- 
eters (Figure  3.1)  were  held  fixed. 

3-U.l  Calculated  Effect  of  Attack 
Angle 

Figure  3.8  shows  the  lateral  acceleration  of  the  projectile  CG  for 
different  initial  values  of  the  attack  angle.  Beyond  the  point  of  nose 
embedment  (0.2  ms),  there  is  a marked  difference  in  acceleration  level 
for  attack  angles  of  1,  3,  and  5 degrees.  The  same  is  true  of  the  bend- 
ing moment  (Figure  3.9). 

An  increase  in  the  initial  attack  angle  clearly  raises  the  magni- 
tudes of  the  lateral  acceleration  and  the  bending  moment.  This  rise  is 
a direct  consequence  of  Equations  2.5,  2.7,  and  2.9  in  Chapter  2.  The 


Appendix  C presents  a comparison  of  lateral  accelerations  obtained 
with  the  WES  and  the  AVCO  formulations  of  the  force  distribution. 
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rate  of  rise,  however,  awaits  experimental  data  for  verification.  If 
future  tests  generally  show  a stronger-than-predicted  increase  in  accel- 
eration and  bending  moment  with  increasing  attack  angle,  then  the  veloc- 
ity terms  in  the  lateral  stress  equations  will  need  to  be  enlarged,  with 

the  a terms  reduced. 
s 

3.4.2  Calculated  Effect  of  Impact 
Velocity 

Aside  from  decreasing  the  lateral  acceleration  and  the  bending  mo- 
ment, a reduction  of  the  impact  velocity  lengthens  the  time  required  for 
penetration,  and  vice  versa  (Figures  3.10  and  3.1l).  From  the  calcula- 
tions, it  also  appears  that  a given  change  in  the  impact  velocity  has  a 
stronger  effect  than  the  same  relative  change  in  the  angle  of  attack. 

In  this  case  a factor  of  2 increase  in  the  impact  velocity  produces 
about  the  same  increase  in  the  lateral  load  as  a factor  of  5 change  in 
the  attack  angle. 

3.4.3  Calculated  Effect  of  Target 
Strength 

Inspection  of  Equations  2.5,  2.7,  2.9,  and  3.3  reveals  that  the 
loads  on  the  projectile  must  exhibit  a target-strength  dependence  that 
falls  somewhere  between  linear  and  square  root.  For  an  impact  velocity 
of  1500  ft/s  and  an  attack  angle  of  3 degrees,  the  lateral  acceleration 
and  the  bending  moment  are  roughly  proportional  to  the  square  root  of 
the  strength,  as  can  be  seen  from  Figures  3.12  and  3.13.  At  much  lower 
velocities,  the  dependence  will  be  more  nearly  linear. 

3.4.4  Calculated  Effect  of  Wake 
Separation 

The  analysis  of  wake  separation  in  Chapter  2 is  insufficient  for 
predicting  separation.  It  is,  however,  adequate  for  assessing  the  con- 
sequences of  a specified  amount  of  separation.  Thus,  by  varying  the 
ratio  6/r  , it  is  possible  to  find  out  the  degree  to  which  separation 
(and  reattachment)  may  influence  the  loads  on  the  projectile. 

The  value  6/r  = 0 represents  no  initial  wake  separation,  but 
parts  of  the  projectile  may  rotate  away  from  the  target,  thereby  losing 
contact.  The  value  6/r  = 1 represents  an  unreasonably  large  amount  of 
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separation  (see  Figure  2.4).  For  rock  and  concrete,  ratios  on  the  order 
of  6/r  = 0.1  seem  reasonable. 

Figures  3.14  and  3.15  show  the  results  of  calculations  made  with 
different  values  of  6/r  . The  variation  of  6/r  ha3  little  influence 
on  the  maximum  lateral  acceleration,  affecting  only  the  pulse  duration. 

The  bending  moment  is  more  strongly  affected,  especially  by  the  larger 
values  of  6/r  . Nevertheless,  the  bending  moment  and  the  lateral  accel- 
eration are  essentially  unchanged  when  6/r  is  increased  from  0.0  to  0.1. 
According  to  the  calculations,  a great  deal  of  separation  (6/r  >0.5)  is 
required  to  produce  a significant  change  in  the  bending  moment. 

3.5  AVCO  REVERSE  BALLISTIC  TESTS 

The  Sandia  RBT's  offer  perhaps  the  best  non-normal  penetration  data 
on  record,  since  the  projectiles  were  fully  instrumented  and  the  targets 
were  large  enough  to  delay  edge  effects  for  a whole  body  length  of  pene- 
tration. However,  the  only  nonzero  angle  of  attack  used  was  3 degrees. 

At  least  one  more  value  is  needed  to  establish  or  verify  a relation  be- 
tween the  projectile  loads  and  the  attack  angle  for  a specified  target. 
AVCO  Corporation  has  conducted  RBT's  in  concrete  for  angles  of  2,  5,  and 
10  degrees. 

3.5-1  Test  Description 

The  projectile  parameters,  target  properties,  and  impact  conditions 

for  the  AVCO  tests  are  given  in  Figures  3-16  and  3.17.  A 5-degree  angle 

was  tested  at  1050  ft/s,  and  all  angles  were  tested  at  1500  ft/s.  All 

5 

the  projectiles  were  instrumented  with  strain  gages,  but  only  the 
1050-ft/s  projectile  was  instrumented  with  accelerometers. 

In  each  of  the  tests  a concrete  cylinder  15  inches  across  and  9 
inches  deep  was  fired  from  a gun,  striking  a stationary  3-inch-diameter 
projectile.  Due  to  the  size  of  the  target  and  the  sonic  speed  (*8000 
ft/s)  in  the  concrete,  the  tests  were  free  of  edge  effects  for  no  more 
than  0.2  ms.^ 

Only  the  strain  gage  S,  shown  in  Figures  3.18  and  3-19»  has  been 
, chosen  for  discussion. 

In  the  Sandia  tests,  edge  effects  were  absent  for  about  1 ms,  or  18 
inches  of  penetration. 
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3.5*2  Accelerometer  Data  for 
1050-ft/s  Test 


Figure  3.20  shows  a comparison  of  the  axial  deceleration  data  with 
the  CET  results  for  the  1050-ft/s  test.  Except  for  the  peaks  in  the 
test  record,  the  calculation  overshoots  the  data  by  at  least  a factor  of 
2 after  the  first  0.4  ms.  Thus,  the  penetration  resistance  of  the  con- 
crete is  either  (l)  less  than  indicated  by  the  strength  or  (2)  affected 
by  the  target  dimensions. 

It  is  likely  that  the  low  penetration  resistance  is  the  result  of  a 
low  ratio  of  target  diameter  to  projectile  diameter  (5  in  this  case). 

In  the  Sandia  tests,  this  ratio  was  30,  and  the  calculated  deceleration 
matched  the  test  data.  The  larger  the  ratio,  the  better  the  target 
simulates  a half-space.  (An  infinite  ratio  is  assumed  in  the 
calculations. ) 

Figure  3.21  shows  a comparison  of  the  calculated  and  measured  lat- 

7 

eral  accelerations  for  the  forward  accelerometer.  The  sharp  spikes  in 
the  calculated  results  sure  due  to  the  integration  scheme  used  in  the 

Q 

code,  and  there  is  so  much  oscillation  in  the  test  data  that  the  com- 
parison is  difficult  to  appraise. 

3.5.3  Strain  Gage  Data  for 
1050-ft/s  Test 

The  measured  strains  for  the  1050-ft/s  test  (Figure  3.22)  are  much 
lower  than  the  calculated  strains  for  the  same  reason  that  the  accelero- 
meter data  (Figure  3.21)  were  low.  Furthermore,  the  onset  of  strain  is 
delayed,  as  in  the  Sandia  tests.  After  the  delay,  the  initial  rise  in 
the  measured  strain  is  comparable  with  the  calculated  strain,  but  the 
rise  rate  drops  after  about  0.2  ms,  probably  due  to  edge  effects. 

3.5.4  Strain  Gage  Data  for 
1500- ft/s  Tests 

Measured  and  calculated  strains  for  the  1500-ft/s  tests  are  shown 
in  Figures  3.23  through  3.25.  Again,  the  onset  of  the  measured  strain 


7 

0 The  rear  lateral  accelerometer  failed  to  operate. 

The  same  is  true  for  all  calculated  results  shown  hereafter. 
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is  delayed,  and  the  observed  strain  levels  are  much  lower  than  calculated. 
The  discrepancy  between  the  data  and  the  CET  analysis  decreases  with  in- 
creasing attack  angle,  but  this  decrease  does  not  necessarily  mean  that 
the  predictions  are  more  accurate  for  higher  attack  angles  in  general. 
Figure  3.26  indicates  that  the  calculated  maximum  strain  is  a linear 

function  of  the  attack  angle,  while  the  measured  maximum  strain  is  non- 
o 

linear.  The  nonlinearity  may  be  due  to  edge  effects,  in  which  case  the 
comparison  made  in  Figure  3.26  is  inappropriate.  Overprediction  of  the 
strain  is  to  be  expected  for  small  targets;  but  if  future  tests  for 
large  targets  exhibit  the  same  type  of  nonlinearity,  then  the  CET  analy- 
sis will  require  some  modification. 


o 

Individual  rapid  oscillations  are  ignored 
strain  levels  from  Figures  3.23  through  3 


in  determining  the  maximum 

25. 
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PROJECTILE 


0.284-SCALE  EARTH  PENETRATOR 

WEIGHT  =9.48  LB 

MAXIMUM  DIAMETER  = 1.9  IN 

MOMENT  OF  INERTIA  = 221.9  LB-IN2  (ABOUT  CG) 

NOSE  SHAPE:  6.00  CRH  OGIVE  (BEVELED  TIP) 

Ep  = 2.8  x 107  PSI 


Figure  3.1  Projectile  parameters,  target  properties,  and  impact 
conditions  for  Sandia  RBT's. 
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NOSE  BASE 

CG 

0.00 

CENTER  OF  GRAVITY 

F 

+0.75 
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-8.38 

REAR  ACCELEROMETER 

S 

-0.89 

STRAIN  GAGE 

CG  B 


-10  -5  0 5 10 

DISTANCE  FROM  CG,  INCHES 


Figure  3-2  Axial  location  of  strain  gage  and  accelerometers 
for  projectile  used  in  Sandia  RBT. 
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AXIAL  DECELERATION,  10 


0 0.2  0.4  0.8  0.8  I 

TIME,  MS 

Figure  3.3  Comparison  of  calculated  and  measured  axial 


decelerations  for  Sandia  RBT  (forward 
accelerometer) . 
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LATERAL  ACCELERATION,  I03  q'S 


Figure  3.5  Comparison  of  calculated  and  measured  axial 
decelerations  for  Sandia  RBT  (rear 
accelerometer) . 
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LATERAL  ACCELERATION,  I03 


Figure  3.6  Comparison  of  calculated  and  measured  lateral 
accelerations  for  Sandia  RBT  (rear 
accelerometer) . 
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AXIAL  STRAIN,  PERCENT 


CET  CALCULATION 

TEST  2 

TEST  3 

NOTE:  IMPACT  VELOCITY  1500  FT/S 
ATTACK  ANGLE  = 3 DEGREES 


0.2  0.4  0.6  0.6  1.0 

TIME,  MS 

Figure  3.7  Comparison  of  calculated  and  measured  axial 
strains  for  Sandia  RBT,  strain  gage  S. 
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LATERAL  ACCELEI 


BENDING  MOMENT,  10s  IN  - LB 


LEGEND 


0 0.2  0.4  0.6  0.6  I 

TIME,  MS 

Figure  3-9  Calculated  bending  moment  for  different  angles  of 
attack,  Sandia  RBT,  section  S. 
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LATERAL  ACCELERATION,  I03 


TIME,  MS 


Figure  3.10  Calculated  CG  lateral  acceleration  for 

different  impact  velocities,  Sandia  RBT. 
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BENDING  MOMENT,  10 


Figure  3.11  Calculated  bending  moment  for  different  impact 
velocities,  Sandia  RBT,  section  S. 
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LATERAL  ACCELEf 


BENDING  MOMENT.  10s  IN  - LB 


LEGEND 


Figure  3.13  Calculated  bending  moment  for  different  values  of 
target  strength,  Sandia  RBT,  section  S. 
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LATERAL  ACCELERATION,  I03  g'S 


Figure  3.1^  Calculated  CG  lateral  acceleration  for 
different  amounts  of  wake  separation, 
Sandia  RBT. 
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PROJECTILE 

0. 47-SCALE  EARTH  PENETRATOR 
WEIGHT  = 38.2  LB 

MAXIMUM  DIAMETER  - 3 IN 
MOMENT  OF  INERTIA  2160  LB-IN2  (ABOUT  CG) 
NOSE  SHAPE:  9.25  CRH  OGIVE  (BEVELED  TIP) 
Ep  = 2.8  x 107  PSI 

1,  = 3.98  IN4 


Figure  3.16  Projectile  parameters,  target  properties,  and  impact 
conditions  for  AVCO  1050-ft/s  RBT. 
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PROJECTILE 


0.47-SCALE  EARTH  PENETRATOR 

WEIGHT  36.6  LB 

MAXIMUM  DIAMETER  3 IN 

MOMENT  OF  INERTIA  1937  LB-IN2  (ABOUT  CG) 

NOSE  SHAPE  9.25  CRH  OGIVE  'BEVELED  TIP) 

Ep  - 2.8  x 107  PSI 

ls  - 3.98  IN4 


DIRECTION 


c i = Cj  - 1935  PSI 


Figure  3.17  Projectile  parameters,  target  properties,  and  impact 
conditions  for  AVCO  1500- ft /s  RBT. 
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SIGNIFICANCE 

B 
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NOSE  BASE 
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-15  -10  -5  0 5 10  15 

OISTANCE  FROM  CG,  INCHES 

Figure  3.18  Axial  location  of  strain  gage  and  accelerometers 
for  projectile  used  in  AVCO  1050-ft/s  RBT. 
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STATION 

LOCATION 

INCHES 

SIGNIFICANCE 
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5.36 

NOSE  BASE 

CG 

0.00 

CENTER  OF  GRAVITY 
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3.28 

STRAIN  GAGE 

Figure  3.19  Axial  location  of  strain  gage  for  projectile  used 
in  AVCO  1500-ft/s  RBT. 


U7 


Figure  3.20  Comparison  of  calculated  and  measured  axial 
decelerations  for  AVCO  1050- ft/s  RBT,  rear 
accelerometer. 


LATERAL  ACCELERATION,  I03  g'S 


Figure  3.21  Comparison  of  calculated  and  measured  lateral 
accelerations  for  AVCO  1050-ft/r.  RBT,  forward 
accelerometer. 
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AXIAL  STRAIN,  PERCENT 


T 


Figure  3.22  Comparison  of  calculated  and  measured  strains 
for  AVCO  1050-ft/s  RBT,  strain  Rage  C. 
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AXIAL  STRAIN,  PERCENT 
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Figure  3.23  Comparison  of  calculated  and  measured  strains 
for  AVCO  2-degree  RBT,  strain  gage  S. 


51 


'"I T 


Figure  3.2k  Comparison  of  calculated  and  measured  strains 
for  AVCO  5-degree  RBT,  strain  gage  S. 
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Figure  3.25  Comparison  of  calculated  and  measured  strains 
for  AVCO  10-degree  RBT,  strain  gage  S. 
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MAXIMUM  AXIAL  STRAIN,  PERCENT 


ATTACK  ANCLE,  DECREES 


Figure  3.26  Comparison  of  calculated  and  measured  strain 
levels  for  AVCO  1500-ft/s  RBT,  strain  gage  S. 
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CHAPTER  U 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  CET  analysis  for  normal  impact  and  penetration  has  been  ex- 
tended to  include  non-normal  impact  for  small  angles  of  attack  (less 
than  10  degrees)  and  hard  targets  (such  as  rock  and  concrete).  The  ex- 
tension was  accomplished  without  introducing  any  new  target  parameters 
(empirical  or  otherwise),  the  primary  assumption  being  that  the  lateral 
and  axial  force  distributions  can  be  specified  separately. 

The  non-normal  CET  analysis  was  installed  in  a modified  version  of 
the  AVCO  DAFL  code,  replacing  AVCO's  definition  of  the  force  distribu- 
tion on  the  projectile.  The  modified  code  was  then  used  to  make  calcula- 
tions for  comparison  with  non-normal  RBT  data  obtained  by  Sandia  Labo- 
ratories and  by  AVCO  Corporation. 

The  calculated  results  reproduced  the  Sandia  data  fairly  well,  ex- 
cept for  a time  delay.1  The  calculations,  however,  overpredicted  the 
AVCO  strain  and  acceleration  data  by  a factor  of  2 to  6,  depending  on 
the  angle  of  attack.  The  AVCO  targets  were  apparently  too  small  to  simu- 
late the  resistance  of  a half-space  for  depths  greater  than  Ho  percent 
of  the  nose  length  or  times  greater  than  0.2  ms.  In  the  AVCO  tests,  the 
ratio  of  target  diameter  to  projectile  diameter  was  5>  in  the  Sandia 
tests,  this  ratio  was  30. 

The  agreement  between  the  calculations  and  the  Sandia  data  is  en- 
couraging, but  more  tests  are  needed  to  check  the  theory  at  angles  of 
attack  other  than  3 degrees.  If  the  trends  observed  in  the  AVCO  tests 
are  borne  out  in  future  tests  with  larger  targets,  then  the  CET  analysis 
may  require  some  modification.  Nevertheless,  the  theory  as  it  stands 
seems  to  be  accurate  enough  for  projectile  design  calculations  in  the 
meantime. 


The  time  delay  is  apparently  caused  by  nonrigid  projectile  response  to 
the  rapidly  varying  applied  loads.  Appendix  B presents  the  results  of 
a dynamic  structural  analysis,  revealing  a similar  internal -response 
time  delay. 


JL. 
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APPENDIX  A 


CAVITY  EXPAND ION  THEORY 

Reference  2 documents  an  analysis  of  quasi-static  spherical  cavity 
expansion  in  an  elastic-plastic  material.1  This  analysis,  known  as  the 
Cavity  Expansion  Theory  (CET),  is  the  basis  of  the  normal -impact  pene- 
tration theory  (Reference  2).  For  non-normal  impact  and  penetration  in  ( 

hard  targets,  the  effects  of  overburden  and  compressibility  are  negligi- 
ble, and  the  basic  assumptions  and  results  of  the  CET  are  those  given 
below. 

Consider  an  infinite  medium  surrounding  a spherical  cavity  that  1 

slowly  expands  from  zero  initial  radius  to  instantaneous  radius  a > 0 . 

2 

Denoting  the  radial  and  circumferential  stresses  by  and  aQ  , re- 

spectively, the  equations  that  describe  the  bilinear  failure  envelope 
HBA  (Figure  A.l)  on  segments  HB  and  BA  are 

Segment  HB:  - oQ  = (o^  + ag)  sin  + 2c^  cos  ^ (A.l) 

Segment  BA:  °r  “ °q  = + °0)  sin  ♦p  + 2cp  cos  ♦p  (A. 2) 

The  quantities  c^  and  <J>^  are,  respectively,  the  values  of  the  cohe- 
sion and  the  internal  friction  angle  for  segment  HB;  c 2 and  <P^  are 
the  corresponding  values  for  segment  BA. 

The  material  surrounding  the  cavity  is  divided  into  separate  elas- 
tic and  plastic  regions,  as  shown  in  Figure  A. 2.  Denoting  the  radial 
position  by  the  symbol  r , the  relation  between  and  in  the 

region  b _<  r _<  h is  given  by  Equation  A.l,  and  the  corresponding  rela-  i 

tion  in  the  region  a <_  r £ b is  given  by  Equation  A. 2.  The  geometry  i 

is  spherically  symmetric,  and  the  material  equation  of  equilibrium  is 


Reference  numbers  refer  to  similarly  numbered  items  in  the  References 
2 at  the  end  of  the  main  text . 

Compressive  stresses  and  strains  are  considered  positive. 


57 


l 


* 


(A. 3) 


Upon  integrating  Equation  A. 3 from  r = a to  r = ® , it  is  found 

that  the  radial  stress  a at  the  cavity  wall  is 

s 


as  = (c2  cot  <J>2  - c1  cot  4'1)(Ii23  “ 2) 

c2  cos  $2(l  + sin  ^ cos  4^(1  + sin  4>2 ) 

+ sin  - sin  4>2 


(A.U) 


where 


12  T3/Y 
X2 


(A. 5) 


U sin  ^ 
1 + sin  <t>. 


(A. 6) 


U sin 
1 + sin  $r 


(A.T) 


^ ■ (I 

(3  - sin  4>1)[c2  cos  $2  - c^l  + cot  sin  $2J 
3 (sin  4»1  - sin  ^2)c1  cot  ^ 


(A. 8) 


(A. 9) 


The  quantity  1^  is  the  material  rigidity  index,  and  e#  is  the  radial 
strain  at  the  onset  of  plastic  behavior  (which  is  the  same  as  the  axial 
strain  to  failure  in  a triaxial  compression  test). 

The  von  Mises  condition  (c^  ■ c2  = c , = ^2  = °)  is  the  simplest 

failure  condition  for  which  Equation  A.U  is  applicable.  In  this  case. 
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T 


the  expression  for  0g  reduces  to 


o = — c(l  + in  Ij 
s 3 1 


(A. 10) 


If  the  material  is  governed  by  the  linear  Mohr-Coulomb  condition 
(c^  = Cg  = c , $ = = <|>  > 0) , then  Equation  A.U  becomes 


o = 
s 


3(1  + sin  $)(lj/3  - 1) 


3 - sin 


c cot  <{> 


(A. 11) 


The  failure  criterion  for  some  materials  may  be  •••  approximated  using 
a bilinear  condition  in  which  > 0 and  $2  « u (von  Mises  limit). 
When  this  situation  arises,  the  expression  for  simplifies  to 


o 

s 


c2(l  + sin  <t1)  - c1  cos  4^ 
sin  4>. 


+ IC2 


In  I 


12 


(A. 12) 


If  I12  i.  1 (when  0 < c.^  < c2  , 0 < $2  < 4^),  then  the  entire 

plastic  region  is  governed  by  segment  HB  of  the  failure  envelope.  In 

this  case,  c is  set  equal  to  c , 4>  is  set  equal  to  <k  , and  o 

1 X s 

is  given  by  Equation  A. 11. 
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Figure  A.l  Bilinear  failure  envelope. 
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Figure  A. 2 


Slowly  expanding  spherical  cavity  in  an 
infinite  elastic-plastic  medium. 


APPENDIX  B 


NONRIGID  PROJECTILE  RESPONSE 

Chapter  3 discusses  briefly  projectile  deformation  under  static 
load,  ignoring  the  complications  introduced  by  a rapidly  changing 
(dynamic)  load.  Although  dynamic  structural-response  analysis  i3  be- 
yond the  scope  of  Chapters  2 and  3,  it  seems  to  explain  in  part  the 
time  delay  in  the  RBT  data. 

T.  Belytschko^"  has  used  the  CET-calculated  external  loads  in  a 
lumped-parameter  computer  code,  calculating  dynamic  projectile  response 
for  the  3-degree  Sandia*  RBT's.  Figures  B.l  and  B.2  show  Belytschko's 
lateral-acceleration  results  compared  with  the  original  rigid-body 
calculations  (Chapter  3).  A time  delay  is  clearly  evident  in  the  non- 
rigid  response  for  the  rear  accelerometer. 

A comparison  of  the  nonrigid  calculations  with  the  test  data 
(Figures  B.3  through  B.5)  shows  that  the  predicted  initial  response  is 
indeed  closer  to  the  measured  initial  response  than  are  the  rigid-body 
calculations.  The  late-time  accelerations,  however,  are  still  out  of 
phase  with  the  data.  The  best  agreement  occurs  for  the  axial  strain 
date  (Figure  B.5). 

These  nonrigid  response  calculations  do  not  fully  reconcile  the 
CET  predictions  with  the  RBT  data.  Nevertheless,  the  appearance  of  a 
time  delay  (due  solely  to  nonrigidity)  in  the  calculations  is  evidence 
that  the  delays  observed  in  Chapter  3 probably  arise  from  the  same 
source. 


Personal  communication  dated  June  1978  from  T.  Belytachko,  Inc., 
Chicago,  111.,  to  R.  Bernard,  Soil  Dynamics  Division,  Soils  and 
Pavements  Laboratory.* 
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LATERAL  ACCELERATION,  10 


Figure  B.2  Comparison  of  rigid  und  nonrigid  projectile-response 
calculations  for  rear  accelerometer,  Sandia  RBT. 
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’ XT- 


Figure  B.'j  Comparison  of  measured  lateral  accelerations 
with  dynamic  structural-response  calculation 
for  forward  accelerometer,  fundi. a RBT. 
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LATERAL  ACCELERATION,  I03 


AXIAL  STRAIN,  PERCENT 


Figure  B. 5 Comparison  of  measured  axial  strains  with  dynamic 
structural-response  calculation  for  strain  p,ap,e 
S,  Sand i a RBT. 


APPENDIX  C 


COMPARISON  OF  CET  AND  AVCO  CALCULATIONS 

At  present,  the  AVCO  DAFL  seems  to  be  the  only  viable  alternative 
to  the  CET  analysis  for  calculating  non-normal  projectile  loads.  Al- 
though the  DAFL  formulation  is  quite  different  from  the  CET  analysis, 
it  is  worthwhile  to  compare  the  predictions  of  the  two  for  the  same 
projectile,  target,  and  impact  conditions. 

D.  Henderson  of  AVCO  has  conducted  an  extensive  parameter  study 
(Reference  7)  with  the  DAFL.  The  baseline  projectile  and  concrete  tar- 
get from  this  study  have  been  chosen  for  comparison  of  DAFL  and  CET 
predictions.  The  projectile  parameters  are  shown  in  Figure  C.l.  The 
nominal  properties  of  the  concrete  are  as  follows:  Y = 5000  psi,  p = 

O /T 

4.66  slug/ft  , and  E = 3 x 10  psi. 

Figure  C.2  compares  the  results  for  maximum  lateral  CG  accelera- 
tion versus  angle  of  attack.  In  these  calculations,  the  CET  analysis 
consistently  overshoots  the  AVCO  results  by  about  50  percent.  The  dis- 
crepancy, however,  is  in  the  absolute  magnitude  of  the  acceleration; 
the  slope  of  the  curves  is  about  the  same  for  the  same  impact  velocity. 

Since  mechanical  properties  are  usually  more  readily  obtainable 
than  are  AVCO's  empirical  target  parameters,  the  CET  analysis  is  the 
easier  method  to  use  for  before-the-fact  calculations  in  untried  tar- 
gets. Nevertheless,  when  used  by  an  experienced  analyst  such  as  Hender- 
son, the  DAFL  has  generally  proved  to  be  capable  of  accurate  predic- 
tions. Both  methods  should  be  used  whenever  possible. 
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APPENDIX  D 


NOTATION 

a Spherical  cavity  radius 

A Surface  area 

Ao  Cross-sectional  area  of  projectile  at  section  S 

b Radial  position  corresponding  to  break  point  on  bilinear 
failure  envelope 

c Material  cohesion 

c^jC^  Value  of  cohesion  describing  bilinear  failure  envelope 

CG  Center  of  gravity 

dA  Differential  surface  area 

dM(z)  Point  distribution  of  mass  along  the  projectile  axis 

E Young's  modulus  for  target 

Ep  Young's  modulus  for  projectile 

f^  Lateral  component  of  internal  force 

fz  Axial  component  of  internal  force 

Fx  Lateral  component  of  external  force 

Fz  Axial  component  of  external  force 

h Radial  location  of  elastic-plastic  interface 

I Mass  moment  of  inertia  about  CG 

CG 

I Area  moment  of  inertia  of  projectile  cross  section  S 

1^  Material  rigidity  index  (=  2E/3Y  for  incompressible  von 
Mises  materials) 

L Projectile  length 

M Projectile  mass 

r Local  aftbody  radius;  radial  position  in  the  CET 


Projectile  radius  at  section  S 


4 


r 

s 

s Axial  distance  from  any  point  on  aftbody  to  base  of  nose 

v Local  lateral  velocity 

x 

v Local  axial  velocity 

z 

V^  Lateral  component  of  CG  velocity 

V Axial  component  of  CG  velocity 

z 

x Lateral  distance  from  CG 

X Horizontal  distance  from  origin  fixed  in  target 

X,Z  Components  of  acceleration  in  the  X-  and  Z-directions , 

respectively 

Y Unconfined  compressive  strength 

z Axial  distance  (forward)  from  CG 

z z-coordinate  of  the  nose  tip 

n 

z0  Lateral  velocity  due  to  pure  rotation 

Z Vertical  distance  from  origin  fixed  in  target 

a Angle  of  attack 

T Bending  moment 

6 Wake-separation  gap  width 

e Axial  strain 

a 

e#  Compressive  radial  strain  at  failure 
n Angle  between  projectile  surface  and  projectile  axis 

0 Angle  of  rotation 

• 

9 Angular  velocity 

0 Angular  acceleration 

5 Variable  of  integration 


o'  ,o"  ,o"  ' 

XXX 


o'  ,ct"  ,o'  ' ' 
z z z 


Vy2 


Normal  stress  due  to  axial  motion 
Radial  stress 

Radial  compressive  stress  at  surface  of  slowly  expanding 
spherical  cavity 

Lateral  component  of  stress 

Contributions  to  lateral  stress  due  to  axial,  lateral, 
and  rotational  motion,  respectively 

Axial  component  of  stress 

Contributions  to  axial  stress  due  to  axial,  lateral, 
and  rotational  motion,  respectively 

Circumferential  stress 

Internal  friction  angle 

Internal  friction  angle  describing  bilinear  failure 
envelope 

Azimuthal  angle,  measured  around  projectile 
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Laboratory  Technical  Publications,  a facsimile  catalog 
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Bernard,  Robert  S 

Non-normal  impact  and  penetration:  analysis  for  hard  targets 
and  small  angles  of  attack  / by  Robert  S.  Bernard,  Daniel  C. 
Creighton.  Vicksburg,  Miss.  : (J.  S.  Waterways  Experiment  Sta- 
tion ; Springfield,  Va . : available  from  National  Technical 
Information  Service,  1978. 

79  p.  : ill.  ; 27  cm.  (Technical  report  - I).  S.  Army  Engineer 
Waterways  Experiment  Station  ; S-78-14} 

Prepared  for  Defense  Nuclear  Agency,  Washington,  D.  ()., 
under  Subtask  Y99QAXSB21 l , Work  Unit  18. 
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